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ABSTRACT 

A portable underwater source-receiver system, 
supporting equipment, and digital computer software have 
been developed to determine near-surface sound speed and 
attenuation at harmonic frequencies from 5 to 95 kHz 
simultaneously, and in-situ from a ship. These data have 
yielded values of bubble densities as a function of bubble 
radius from 40 to 1000 microns, depth to 15.2 meters, time 
of day, varying oceanographic conditions and environmental 
parameters for two days in November in Monterey Bay. 
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I. INTRODUCTION 



A. BACKGROUND 

When a sound wave propagates through the ocean, its 
intensity is weakened with increasing distance. There are 
two mechanisms which adequately explain this weakening: 
divergence and absorption. If we are given an infinite 
medium (i.e., no boundaries) that is also homogeneous and 
free from scatterers, this weakening is predictable. 
Spherical spreading changes the sound intensity in propor- 
tion to the inverse square of the range. Absorption 
weakens the sound as a function of frequency. 

In the real ocean we can approach a test situation 
where boundary effects can be ignored. However, the nature 
of the sea makes it impossible to assume a homogeneous 
medium. It has been found through observations of sound 
waves in the ocean that both their amplitudes and phases 
fluctuate. This is to be expected since the ocean is, in 
fact, dynamic and is moved constantly by tide and current. 
As this transport of water past a fixed point occurs, the 
inhomogeneities also move, thereby altering the acoustic 
conditions along the sound path. In- situ measurements of 
the ocean microstructure are difficult to obtain by tradi- 
tional methods. 

Acoustic probing, since it does not alter the medium is 
an excellent research tool for learning about the structure 
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of the sea. This probing can be used to determine the 
presence of scatterers by observing their effect on both 
sound attenuation and sound speed as a function of frequency. 

The presence of bubbles in the near-surface ocean causes 
additional attenuation and makes the phase speed dispersive. 
The near-surface frequency dependence of sound propagation 
is especially important with respect to naval sonars which 
often transmit their sound in a surface duct mode. This 
transmission could be subjected to unexpected errors if a 
bubble-free model of computing sound speed were used in an 
area containing a large population of resonant bubbles. 

If a close correlation can be found between bubble popula- 
tion density and weather conditions or environmental param- 
eters, this would enable the Navy to publish sound speed 
and attenuation corrections to sonar performance tables 
which could be applied aboard the individual ship according 
to the existing conditions. 

B. OBJECTIVE 

The objective of this research was to take in- situ 
acoustic measurements over a fixed range of approximately 
one meter in the upper fifteen meters of the ocean in order 
to observe and record sound wave amplitudes and sound 
speeds as a function of frequency. These data were then 
reduced to compute bubble populations by size and depth 
and to compute a differential value of sound speed as a 
function of frequency in the upper reaches of the sea. 
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In order to make the observations, a constant amplitude 
sawtooth sound wave, rich in harmonic content, was trans- 
mitted horizontally through the sea past two acoustic 
receivers that were spatially fixed in the far field of the 
transducer. Signals received by these hydrophones were 
digitized and processed by a high-speed, high-resolution 
computer. Amplitudes recorded at various depths and loca- 
tions, at various times, and under various oceanographic 
and meteorological conditions were spectrally analyzed, 
and were then related to bubble population of the medium. 
The computed phase angle of each harmonic permitted the 
calculation of the sound speed at that frequency. The 
computed propagation loss, compared to the bubble-free 
medium yielded the attenuation at that frequency. 

C. BUBBLE ABSORPTION AND SCATTER THEORY 

Bubbles have a drastic effect on amplitude attenuation 
if they are of such a size that they are at or near the 
resonance frequency of a transmitted component. If this 
situation exists, the bubbles resonate, thus absorbing 
and scattering that particular frequency component of the 
incident sound wave, and causing changes in the sound 
speed. It is this property of bubbles that was utilized 
in this research project to determine average bubble popu- 
lations over the sound path. 

Determination of bubble population by acoustic methods 
is shown to be theoretically practical by considering the 
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scattering cross section, , and the absorption cross 
section, a , of a bubble. Computation shows that, at 
resonance, the scattering cross section of the most 
prominent bubbles is 180 to 900 times greater than the 
geometrical cross section and the absorption cross section 
is many times again larger than the scattering cross sec- 
tion [Medwin, 1970] . Because of this inordinately large 
effective bubble cross section at resonance in conjunction 
with an extremely narrow resonance peak, it is then quite 
easy to detect resonant bubbles even in the presence of 
non-resonant bubbles, fish, and particulate matter. 

The size distribution of bubbles has been investigated 
by numerous research scientists. Most notably, efforts 
by Glotov [1962] and by Blanchard and Woodcock [1957] imply 
that the majority of bubbles encountered are below 200 
microns radius with a population peak in the vicinity of 
60 microns. Larger bubbles tend to rise rapidly to the 
surface and disappear due to buoyancy. Smaller bubbles 
tend to shrink and disappear. This information indicates 
that a frequency range from 10 kHz to 100 kHz should effec- 
tively define the major bubble population. For this 
experiment, a range from 5 kHz to 100 kHz was chosen. 

For present purposes, let us assume bubbles of a uniform 
size are present in a sound field. Further assume no 
bubbles of other sizes present, and that the number is 
small enough that interaction may be discounted. The 
attenuation due to the presence of these bubbles is then 
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the additive contribution caused by the extinction effect 
of each individual bubble. The extinction effect is the 
combined effect of absorption and scattering, that is, the 
extinction cross section, cr , is the sum of the scattering 
cross section, o , and the absorption cross section, cr 

S cl 

[Class Notes, Medwin, 1964]. 



where 



cr = cr + cr , 
e s a 



( 1 ) 



cr 



Power scattered over all angles 
Incident Acoustic Intensity 



47rR o 

((to /wW-i)^+6 : 
o 



( 2 ) 



cr = 



Power absorbed by bubble 
Incident Acoustic Intensity 



o (6/K R -1) 
s oo 



(3) 



Therefore, after consolidation 
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( (to Q / to) z - 1 ) Z +6 Z ^- 


o 

to 

o 



(4) 



where 

Rq = Resonant bubble radius, 

to = Angular frequency of incident sound, 

Wq = Angular frequency of bubble pulsation at resonance, 

6 = Total bubble damping constant, 

K n = ^ = Sound propagation constant at bubble 
u \ Q 

resonance frequency. 
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If the intensity incident on the sound field is Iq, then 

the amount of power dissipated by an individual bubble is 

la . If we let N(R) equal the number of resonant bubbles 
0 © 

present in a unit volume, then the change in intensity over 
a distance dx is 

. "OpN (R) x , /c\ 

dl = -I a N (R) e e dx (5) 



o e 



after integration 



I = l 0 e 



-a e N (R) x 



( 6 ) 



Expressing the change in intensity level, or attenuation, 
in decibels 

I 



AI . L . = 10 Log 



10 I 



= -10a N(R)xLog in e 
e i U 



= 4.34 a N (R) x 
e 



(7) 



Finally, expressing this attenuation relative to a unit 
distance , 

AI.L. = 4.34 a N (R) (8) 

e 

x 

In-situ measurements of bubbles at sea, of course, involve 
bubbles of various sizes. It becomes necessary then, to 
determine the degree of response of near-resonant bubbles 
and the range of bubble sizes that will cause sound of a 
unique frequency to attenuate significantly over a known 
distance. From previous studies, it has been found that 
appreciable scattering and absorption will occur only for 
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frequencies relatively close to the resonance frequency. 

In fact, it has been found that if the frequency range, Af, 



is equal to fg ±0.045 times the resonant frequency, then 
the extinction cross section is only one half its resonance 
value [Class Notes, Medwin, 1964]. Thus, it follows that 
only bubbles of size quite close to resonant size will 
contribute to the extinction cross section. 

To calculate this total extinction cross section, it 
will be necessary to integrate over a range of bubble 
radii. Theoretically, the integral will be over the entire 
range from zero to infinity. Practically, however, it has 
been shown that only a small range of sizes need be con- 
sidered in close proximity to the resonant bubble size. 
Thus, our total extinction cross section, o^, reduces to 



S E " 



R +dR 

0 O n (R) dR = 
e 



rR+dR 4TTR 2 (6/6r)n(R)dR 

T(5 /w) 2 -i) 2 + s 2 (9) 

R 0 



where 

n(R)dR = Total number of bubbles per unit volume with 

radii between Rq and R^+dR, 

<S r = Resonant bubble damping constant, 

_ ? 

= KqRq (1-36x10 for a clean air bubble in 
water at one atmosphere. 

In view of the small range of radii considered for each 
resonance frequency, it is valid to assume the quantities 
n(R) , 6, and <5^ are constant over the range and, therefore. 
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treat them external to the integral. This yields 



S_ = 4tt (6/ 6 )n (R) 
sl r 



/•R +dR 
o 

R 

J 0 



2 

R dR 

( (&> o /w) 2 + 5 2 



( 10 ) 



Next, it is necessary to express all variables in the 
integral in terms of a common variable. 

Let 6? = Uq/o) -1 = R/Rq -1. 

Then d(K = dR/R Q . 

Substituting 



s 



E 



= 4TT(6/6 r )n(R) 



'+(R 



(CR + 1) (Rn)) 2 Rnd(R 

(((J} + 1) 2 -1) 2 +6 2 



And, after consolidating terms 



S 



E 



4ttR 3 (6/ 6 ) n (R) 
o r 



'+(R 

.-(R 



q+R) 2 d(R 

„ 9 2 2 

((l-HR)--i) +6 



(ID 



Now, simplifying by dropping higher order terms in(^ , we 
arrive at 



S 4ttR (6/6 n (R ) 
E o r 






d(R 



2 „ 2 



-R 4^ +6 



And, after integrating 



S E ^ 4irR Q (6 / 6 r ) n (R ) 



1/ 6 Tan 1 2(55/6 



= 4ttR 3 n(R)/5 Tan 1 2(R / 6 
o r 



( 12 ) 



17 



The maximum extinction cross section will occur when 
Tan" 1 2fc/<5 = tr/2. If we choose reasonable, representative 
values for 6^ and 6 such as =0.1 and 6 = 0.1, Tan 26^/6 
= 0.71(iv/2). Thus, bubbles within 10% of the resonant 
bubble size will cause 71% of the total extinction cross 
section. 

For our purposes, and in order to predict the maximum 
extinction cross section with maximum attenuation, we will 
evaluate in the limit with Tan 1 2G(/ <5 = tt/2. Thus, equation 
(12) reduces to 

S £ = 2ir 2 R o 3 n (R) / <S r (13) 

Returning to equation(8), the expression equating attenuation 
due to single size bubbles, we may now substitute our total 
extinction cross section, S £ , for a e N(R), the extinction 
cross section for single radius bubbles. This substitution 
yields an expression for attenuation due to all bubble sizes, 
a = 4.34 S E 

= 4.34 2tr 2 R o 3 n(R) /6 r (14) 

For a meaningful method of expressing the attenuation m 
terms of the volume of air in the bubbles, we can define 
U(R)dR as the volume of air per unit volume of water for 
bubbles with radii between R Q and R Q +dR. Assuming spherical 
bubbles, this quantity can be expressed 

U(R)dR = 4/3 TrR Q 3 n(R)dR (15) 
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3 

Cancelling and solving for Rq n(R) yields 

R 3 n(R) = 3/4 ttU(R) (16) 

o 

Substituting this expression into equation (14) and assuming 

- ? 

a clean air bubble with <$ = 1.36x10 

r 

a = 1.50xl0 3 U(R) dB/meter (17) 

The next consideration of major importance is to calculate 
the size of a resonant bubble. The resonant relationship 
is [Eckart, 1945] 

f o ■ 1/2,1 R o /3lfP o /p o (18) 

where y = Effective ratio of specific heats of bubble gas, 

Pq = Density of surrounding water, 

Pq = Ambient pressure at bubble depth, 

t = Measure of surface tension. 

t deviates only slightly from unity and then only for very 

small bubble radii, and thus can be discounted [Devin, 1959] 

y for clean air bubbles falls in the range 1.16 < y £ 1.40 

for bubbles larger than 20 microns. For prediction purposes 

the value 1.40 will be used. This will cause an error no 

greater than 10% in resonant bubble radius. Pq for sea 

water at 20° C and standard atmospheric pressure is 
3 

1026 Kg/M . Finally, ambient pressure at bubble depth 
may be calculated from the linear expression 
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p 






depth (meters) 
1+ 10.337 meters 



NT 



(19) 



1.0x10 



5 



Thus, equation (18) reduces to 



1 .02x10 



-2 



R 



f 



/“P 



( 20 ) 



o 



o meters 



o 



Use of this equation yields TABLE 1, a compilation of 
resonant bubble radii as a function of frequency and depth, 
from the surface to 16 meters. 

In order to demonstrate this theoretical approach, a 
sample computation will best demonstrate the application of 
the resonant bubble theory. 

Given: 60 kHz sound source 



Matched receivers at 5 meter depth 
Unbounded salt water medium 
One meter spatial path between receivers 
Data yields excess attenuation from receiver 
one to receiver two of 1.00 dB relative 
to 15.2 meter bubble-free calibration depth 



P Q at 5 meter depth 



P 



o 



5 ( 5 NT 

1,0x10 1+ io.337i M 7 



NT 



P 



0 




Rq at 5 meter depth for 60 kHz signal 



R - 1/2* f 0 3 Y P 0 /p 0 
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= 1 / 2 TT 



|(60xl0 3 ) 



3(1.4) ( 1 . 5xl0 5 / 1026) 



h 



R =66 microns 

o 



a = l.SxlO 3 U (R) 



U (R) = 1.0/1.5xl0 3 



U (R) = 6.7x10 



-1 

m 



Finally, solving equation (16) for the number of resonant 
bubbles in a 1 micron band 

n(R)dR = 3/4tt U(R)dR (1/R q 3 ) 

n (R) dR = 3/4 ( 6 . 7 x 1 0" 4 ) ( 1 . Ox 1 0 _6 ) ( 1 / 6 6x 1 0 _6 ) 3 



n(R)dR = 556 resonant bubbles per cubic 



meter of radius 66 ±0.5 microns 
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35.1 
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38.0 


39.2 


40.5 



RESONANT BUBBLE RADII AS A FUNCTION 
OF FREQUENCY AND DEPTH 



Assumptions : Clean Bubble Surface 

Y = 1.40 

No Surface Tension Effects 
No Heat Conductivity Effects 



TABLE 1 
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D. SOUND SPEED DISPERSION THEORY 

The sound speed dispersion theory was described in 
Rautmann [1971] and Fitzgerald [1972] , portions of which 
are included here for information. 

The following equation relates bubble size to the speed 
of sound and indicates the speed dispersion for a bubble 
population of unique radius R 
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( 21 ) 



C Q = nondispers ive sound speed, 

o) = sound frequency, 

o) Q = resonant frequency of bubble of radius R , 

U(R ) = fractional air-to-water volume of bubbles having 

° radius R , 
o * 

6 = damping constant of bubbles of radius R . 

This relation predicts that in a region of bubbles with 
a unique radius and a resonance frequency io , the sound 
speed will increase if an incident sound wave has a fre- 
quency greater than o) , and will decrease if the incident 
frequency is less than o) . Figure 1 shows typical sound 
speed changes in an area of unique bubble radii. 
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FIG. I EFFECT ON SOUND SPEED OF A UNIQUE RADIUS BUBBLE POPULATION 



For frequencies far below the resonance frequencies of 
bubbles present in a mixture, w << co q , equation 21 reduces 
to 



c = c 

o 



1 



3U 

2(kR ) 2 
o 



( 22 ) 



A bubble population including several different bubble 
sizes would have the same equation predicting the speed of 
sound, except the fractional volumes and resonant frequency 
effects would have to be summed over all bubble radii 
present, causing a combination of several speed variation 
curves such as Figure 1, with the amplitudes of each pro- 
portional to the density of that particular sized bubble's 
population. 

If a region of the ocean contained an unusually large 
population of bubbles near one resonant bubble size in 
comparison to the population of bubbles within ±20% of 
that size, a curve such as Figure 1 would be expected. 

Since bubbles have little effect on frequencies more than 
20% removed from their resonance frequency, more than one 
such dispersion area could exist on the same curve centered 
on different frequencies. This would represent more than 
one unusually large population of bubbles centered around 
those resonant frequencies. 

The computational method for determination of the speed 
of sound using the phase difference was the same as that 
used by Rautmann [1971] once the phase angle of each 
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component harmonic was given by either the Fast Fourier 
Transform or by using the phase meter. This method is 
repeated here. 

Each harmonic frequency received by the first 
hydrophone can be written as: 

y A = a sin (wt) (23) 

and the signal received by the second hydrophone is: 



y fi = b sin (wt-kx) 



(24) 



where 

k = the wave number = ^ , 

f = the frequency of the applied signal, 

C = the wave phase speed, 

x = the separation between the hydrophones. 



Since we compare the two signals with respect to phase 
by using the digital phase meter or the FFT output we have 
to find values for which kx equals the total phase difference 
between the two signals due to an integral number of wave 
lengths and a fraction of the wave length. This is written 
as : 



kx “ (2Tr) n + 

kx = ^ x =(2Tr)n + 0^_ 

r - MX 

(2Tr)n + 0 r 



n = an integer 

0 r = fractional wave length 
in radians 

(25) 
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The final result is: 



c 



f x 



n + 



0 

360 



where 0 is in degrees. 



( 26 ) 



The integer n can be pre-estimated by calculating the number 
of wavelengths using the speed of sound given by the 
velocimeter . 



II. THE EXPERIMENT 



The theory of resonant bubble attenuation has been 
supported by previous theses. However, past research has 
been limited to analog time averaging of single frequencies 
with a sensor system suspended from a stationary platform. 
This research project has attempted to devise a system that 
can be effectively employed from a moving platform, in this 
instance, the Research Vessel ACANIA. Further, a high-speed 
digital processing system has been added to allow simulta- 
neous real time digital analysis of 20 harmonics of a 
sawtooth waveform. 

A brief and simple description of the experiment sets 
the scene for detailed discussion of each facet of the 
project. An underwater test assembly was suspended from 
the Research Vessel ACANIA' s fantail hydraulic crane. The 
test rack was lowered to a depth of 15.2 meters, a depth 
determined from previous research to be bubble - free 
[Medwin e_t aJ . , 1975]. Test rack depth was determined by 
pulsing through an upward looking transducer. This pulse 
was reflected from the surface and received on either of 
the hydrophone receivers mounted on the test rack. By 
triggering an oscilloscope on the transmitted pulse, the 
time required to reach the hydrophone was displayed on the 
CRT. This time was then halved to give one-way travel 
time and converted to depth by multiplication with the 
nominal 1500 meters/sec speed of sound through sea water. 
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When signal processing was accomplished by digital 
methods, a continuous sawtooth waveform, rich in harmonic 
content, at a fundamental frequency of 5 kHz, was trans- 
mitted into the water. This sound wave was directed past 
two fixed hydrophones approximately one meter apart. 

Figure 2 shows the CRT display of this transmitted signal 
along with the received signals after amplification and 
filtration. The amplified and filtered outputs of these 
hydrophones were sampled with an analog/digital converter 
and fed to a high-speed digital computer. This signal 
processor computed and displayed the attenuation as well 
as the fractional phase differential between the two 
hydrophones . 

When analog methods were used to obtain data, pure 
sinusoidal waveforms were transmitted into the ocean. The 
returning signals after amplification and filtration, were 
passed to a vacuum-tube voltmeter for amplitude measurement 
and to a phasemeter for relative phase measurement. 

The 15.2 meter depth run was, then, a calibration run. 

As the test rack was raised closer to the surface and 
additional runs were made, the excess attenuation and phase 
shift above that measured at 15.2 meters could be attributed 
to resonant bubbles present in the sound field. Application 
of the theory discussed in Section I-C of this thesis 
permits computation of bubble populations as a function of 
size and depth. Further application of theory contained in 
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Section I-E enables caluclation of the change in speed and 
sound directly attributable to the presence of these 
resonant bubbles. 
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Ill . UNDERWATER EQUIPMENT 



A. SOURCE TRANSDUCER 

The source was the USRD F-33 directional transducer 
designed for operation from 1 kHz to 150 kHz. It has two 
arrays, an inner array of lead zirconate-lead titanate 
ceramic elements measuring 1.5" by 2.0", and an outer 
array of barium titanate ceramic elements measuring 8.0" 
by 8.5". These dimensions result in a far field, a /A, 
of 0.76 meters from the transducer face at 100 kHz. The 
transducer is designed to take a maximum voltage input of 
200 volts (rms) . The typical transmitting response ranges 
from a pressure of 17 dB re 1 ubar/volt at 5 kHz to 48 dB 
re 1 ubar/volt at 100 kHz. These values are computed for 
a distance of one meter and assume 100 feet of attached 
cable . 

B. RECEIVING ELEMENTS 

The hydrophones chosen were Atlantic Research Corp. 
Model LC-lOs. These hydrophones exhibit a flat frequency 
response over the frequency range 1 kHz to 100 kHz with a 
voltage sensitivity in the vicinity of -108 dB re 1 volt/ 
ubar. The sensing element is of lead zirconate titanate 
encased in a bonded neoprene sheath 1.13 inches in length 
and 0.38 inches in diameter. The clement is then connected 
by six feet of low noise cable to a Model 2010 NUS 30 dB 
preamplifier encased in a sealed, waterproof can. This 
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preamplifier has a gain accuracy ±0.2 dB at 1 kHz, is usable 
to 200 kHz, and is powered by 12 volts DC at 4.5 milliamps. 

When transmitting in the pulse mode to check for frame 
reflections, a definite pulse was seen on the receiving 
equipment with no phase lag from the transmitted signal. 

This proved to be an electromagnetic field that was 
generated simultaneously with the sound generation. 

Acoustic signal to electromagnetic noise ratios were found 
at all frequencies of interest, and showed this noise to 
be more than 30 dB down from acoustic signal levels. 
Geometric phasor diagrams at one hydrophone showed this to 
be the equivalent to an error of ±3.2% in attenuation when 
the signals were in phase or 180° out of phase. The phase 
angle effect was most pronounced when the phasors were 90° 
out of phase, causing an angular error of ±1.8°. If the 
signal, in traveling to the second hydrophone, should have 
a similar relation but in the other direction, the combined 
effects would be 6.4% attenuation error and 3.6° phase 
error. The phase angle changes caused by variations in the 
speed of sound during a normal experimental run were 
sufficiently small that most of this error would have 
canceled out in the subtraction of the calibration run 
data, making the expected errors acceptably small. 

When transmitting a sawtooth continuous waveform, 
however, a very significant electromagnetic signal was 
transmitted at each sharp change in voltage direction. 

Since it was impossible to accurately determine which 
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harmonic frequencies were contained in this signal, no 
meaningful signal to electromagnetic noise ratios could 
be calculated, and it became necessary to eliminate the 
radiation altogether, or reduce it to 30 dB below the 
smallest harmonic component received. The received 
electromagnetic signal was finally reduced to undetectable 
magnitudes by completely shielding the LC-10 hydrophone 
and transmitting cable from the sensing element to the 
preamplifiers. This shielding was done with copper mesh 
braided shielding which was grounded to the preamplifier 
casing. This shielding was then insulated from the water 
with neoprene and tape, and the sensing element was insu- 
lated with a coating of epoxy. The sensitivity loss 
caused by this shielding was approximately 3 dB , but since 
it was constant under all sampling conditions, it also 
canceled out in subtraction of the calibration run. A 
previous design, using neoprene for covering the sensing 
element showed evidence of bubble entrapment. This was 
avoided by using epoxy which had been heated and degassed 
by pumping prior to application to the sensing element. 

C . SUPPORT FRAME 

The initial design of the underwater test rack was 
taken from previous thesis work with minor modifications. 
Figure 3 shows, in schematic form, this test assembly with 
all modifications made during test and development. The 
frame was constructed of 1-1/4" hollow steel pipe. In 
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FIG. -3 . UNDERWATER TEST RACK AND ASSEMBLY MODEL I 



order to minimize frame reflectivity, the pipe frame was 
sleeved with heat-shrunk rubber and silicon rubber wedges 
were built up inboard on the horizontal arms of the frame. 
The source transducer was mounted midway on the vertical 
frame member. Heavy duty nylon shotline supports drawn 
between rubber shock mounts were affixed parallel to the 
vertical frame member, designed to support the passive 
hydrophone elements. Sealed cans containing preamplifiers 
for these hydrophones were mounted outboard on the upper 
horizontal frame member far back of the sound field. An 
upward looking transducer was mounted behind these pre- 
amplifiers to be used to determine depth. On the underside 
of the upper arm, a current probe was mounted for remote 
monitoring of the water flow past the test assembly. 
Finally, the reference velocimeter was mounted in the 
vertical position and suspended from the upper horizontal 
extension behind the source transducer. 

Fully assembled, with all equipment mounted in place, 
this test assembly weighed approximately 150 pounds. 

Figure 4 shows this test assembly with all supporting 
equipment in place. 

From previous research work, it was known that to 
obtain acceptable accuracy both in amplitude and phase, 
it was necessary to reduce frame reflections to a level 
30 dB below the direct path amplitude. The test assembly, 
less velocimeter, was lowered into a fresh water anechoic 
laboratory tank for determination of reflection magnitudes. 
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FIG. 4 RACK MODEL I WITH SUPPORTING EQUIPMENT 
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This testing was done using a gated pulse. By carefully 
measuring both direct and specular reflection path lengths, 
and measuring the relative amplitudes of the direct and 
reflected received signals on an oscilloscope, it was 
possible to identify each of the path contributions uniquely 
based upon travel time along the path. As was true with 
electromagnetic radiation contributions, this -30 dB high 
limit on reflection amplitudes resulted in a maximum attenu- 
ation error of 6.4 % and a maximum phase differential error 
of 3.6°. Again, however, the small temperature differences 
between the 50 foot calibration depth and the upper ocean 
depths of interest ensured the speed of sound and there- 
fore the acoustic wavelength would change by such a small 
amount that most of the error induced by reflections would 
cancel out in subtraction. The goal of 30 dB direct 
signal to reflected signal ratio was attained at all 
frequencies of interest on this original test assembly. 

First attempts to use this test assembly configuration 
at sea were totally unsuccessful. Motion of the apparatus 
through the sea imparted by the floating platform caused 
instability problems. Transient noise signals injected 
by snapping of the primary rack support cable totally 
obliterated the received signal. Apparent relative move- 
ment and vibration of the hydrophone support lines added 
to this instability both in amplitude and phase. Attempts 
to minimize these problems by mounting flexible wire 
strong-backs to the assembly for frame stability and 
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adding a fine horizontal support line between the hydrophone 
elements to reduce relative movement met with only partial 
success. The test assembly remained unsatisfactory with 
each of these modifications. Due to the nature of the 
difficulties encountered, it did not appear feasible to 
further modify this test rack design with any expected 
degree of success. The obvious design requirements indi- 
cated the need for a more stable test vehicle. 

The problem of support frame stability was confronted 
first. It appeared most desirable to increase the rigidity 
of the design. In making a more rigid and compact frame, 
reflectivity loomed as the largest obstacle. Not only 
would the added surface area of the reflector and its 
closer proximity to the sound field cause difficulties, but 
the pulsation method of determining reflective contribu- 
tions would be unusable because of the small differential 
in path length and hence signal travel time. It was 
decided that this problem could be circumvented by making 
use of the 15.2 meter bubble-free calibration run. It 
was hoped that in computing the excess attenuation at the 
upper ocean depths by subtracting out the reference attenu- 
ation at 15.2 meters, reflective contributions could be 
lost in cancellation. 

Figure 5 shows schematically the structure of test 
assembly Model II. All support equipment previously used 
on rack Model I was adapted for use on assembly Model II. 
Longitudinal construction was 2-3/8" PVC pipe. Mounting 
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brackets for the source transducer and the reference 
velocimeter were made of 3/8" aluminum. Mounting brackets 
for hydrophone element supports were of 1" lucite. These 
support brackets further provided frame rigidity. Total 
weight of the fully-rigged test assembly was approximately 
150 pounds. 

In spite of the superior performance of assembly 
Model II as regards signal stability and frame rigidity, 
analysis of at-sea data indicated that reflective inter- 
ference from Model II was indeed excessive. Assembly 
Model I was employed for the remainder of the research 
work. Nylon fish line under slight tension (dashed lines 
in Figure 3) was added between the hydrophone support lines 
in an attempt to minimize relative movement. Additional 
fish line was used between the frame and the support lines 
for vibration damping. Cushioning was added between the 
preamplifier canisters and the support frame to eliminate 
vibration transmission. These devices helped somewhat; 
however, the assembly Model I still was usable only in 
relatively mild sea and weather conditions, and there was 
continuing evidence of incorrect phase measurements at the 
lower harmonics of the fundamental sawtooth frequency. 
Finally, stuffing the tubes with steel wool for vibration 
damping made the assembly usable in moderately heavy sea 
conditions . 
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D. SHIP MOTION ISOLATION EQUIPMENT 

The next problem was to minimize the transient signals 
imparted to the rack by the main support cable. The solu- 
tion was to minimize the submerged weight and to support 
the test assembly from a non-linear elastic support line 
creating a simple harmonic oscillator with a period which 
would be relatively uninfluenced by the motion imparted 
to the ship by passing waves. 

A 20-foot section of one-quarter inch cotton- impregnated 
elastic cord was obtained, which had a low constant of 
elasticity (s = force/distance) of 0.029 lbs/inch when 
the applied force was approximately five pounds. This 
corresponds to a resonant frequency ( w n = /s/m) of 0.273 
rad/sec and a period of 23 seconds assuming an assembly 
mass of 150/g slugs. For all ocean swell waves of period 
about 11 seconds, the system was relatively unaffected by 
passing waves. 

To achieve the five pound weight necessary to arrive 
at this low resonance frequency, a method of flotation 
was necessary. This flotation device was required to 
achieve a negative buoyancy of five pounds when the test 
assembly was submerged at various depths. This device had 
the further advantage of increasing the mass of the assembly 
thereby increasing the natural period of oscillation as well 
as providing a damping effect that further isolated the 
assembly from surface motion. 



44 



A cylindrical container, open at one end, was 
constructed of aluminum and was filled with varying 
amounts of styrofoam in an attempt to obtain the desired 
buoyancy. The compressability of the styrofoam made this 
an unusable design for the multitude of depths to be 
investigated. 

In order to provide a constant buoyancy at varying 
depths, a second device was constructed using an air 
bladder which was connected to an air supply on board the 
ship by a flexible air hose. This bladder, which consisted 
of an inner tube within a tire for protection, but open 
to the water, was connected to another fully inflated tire 
mounted on a rim, which acted as a fixed size buoyant 
device as well as a damper. The air pressure in the fully 
inflated tire was kept in excess of the depth pressure 
anticipated to preclude its compression in size. The 
pressure in the bladder was always apprxoimately that of 
the surrounding medium, and increased with depth. Since 
a constant displacement was needed, air had to be added 
with increased depth. 

Figure 6 shows a sketch of the flotation device in 
position above rack Model II. 

E. DEPTH DETERMINATION 

The transducer selected for depth determination was 
the USRD Type F-41 directional transducer designed for use 
from 15 kHz to 150 kHz. The transducer is composed of an 
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array of lead zirconate-lead titanate ceramic elements 
measuring 1.5" by 2.0" that will operate at input voltage 
up to 200 volts (rms). A typical transmitting response 
shows a pressure of 60 dB re 1 ubar/volt at 150 kHz. This 
is computed for a distance of one meter from the transducer 
face and assumes a 100 foot cable. 

F. CURRENT FLOW METER 

Although not necessary as a data input to final 
numerical analysis, a Marsh-McBirney Model 711 Electromag- 
netic Water Current Meter was mounted on the test rack to 
provide oceanographic background for the experiment. This 
device measures water flow in two directions orthogonal 
to the longitudinal axis of the instrument probe. Maximum 
measurable flow is ±10 ft/sec. Accuracy is to within ± 2 % 
of the reading. 

G. REFERENCE VELOCIMETER 

Finally, the velocimeter used to determine the 
reference speed of sound for investigation of the change in 
speed of sound due to resonant bubbles was the NUS Model 
1630-002 SV/T/D sensor. This unit determines sound velocity, 
temperature, and depth. The per cent of error in depth 
readings for the shallow waters of interest rendered this 
section of the instrument useless. The sound velocity 
module gave accuracy to ±0.15 mcters/sec. The temperature 
module was accurate to 0.1° C. 
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IV. SHIPBOARD/ LABORATORY EQUIPMENT 



Several electronic devices were used in various 
configurations depending on the information to be obtained. 
All equipment used will be described, followed by a descrip- 
tion of the set-ups for each experimental arrangement. 

A. SIGNAL GENERATING EQUIPMENT 

To generate a monofrequency signal for the upward 
looking sonar, and for precise triggering of the sawtooth 
function generator, a General Radio Model 1162A Frequency 
synthesizer was used. This provided a highly stable sine 
wave with a maximum output voltage of 2 volts (rms) . 
Frequency synthesizing was accomplished with a quartz 
crystal oscillator. The frequency was digitally selectable 
to 1 MHz with accuracy to 0.01 Hz. 

To generate a monofrequency pulse for determining 
depth, or measuring reflected sound and electromagnetic 
radiation, the synthesizer wave form was gated through a 
General Radio 1396-A tone-Burst Generator. This generator 
permitted the signal to pass or not pass depending on the 
selected number of cycles to be counted. 

To generate a sawtooth waveform for the transducer, 
the synthesizer was used to trigger a Wavetek 144 Sweep 
Generator. The experiment consistently used a 5 kHz tri- 
angular waveform with a 19:1 asymmetry as an approximation 
of a sawtooth wave. 
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After generation, the output voltage was amplified by 
a Hewlett-Packard 467A Power Amplifier to a maximum of 
±20 volts (peak) before being passed to the appropriate 
transducer . 

A Fluke 8000A digital multimeter was used to insure 
precise amplitudes of transmitted voltages. 

B. SIGNAL PROCESSING EQUIPMENT 

All signals were received by the LC-10 hydrophones and 
amplified* by the 30 dB NUS preamplifiers as described in 
Section III-B. They were then transmitted through 150 feet 
of shielded cable to the ship where they were evaluated. 

The received signals were immediately filtered to reduce 
low frequency motion- induced and ambient noise, and to pre- 
vent high frequency aliasing as discussed in Section V-A. 

The filtration was done using Krohn-Hite 3550(R) filters 
with a bandpass adjustable to any bandwidth, and a range 
between 2 Hz and 200 kHz, with a roll-off at both the high- 
pass and low-pass frequencies of 24 dB/octave. These 
generated less than 200 microvolts of noise and could 
handle a signal up to ±7 volts (peak) . 

The filtered signals were further amplified by a pair 
of Princeton Applied Research (PAR) 113 Low Noise Preampli- 
fiers which could amplify signals with frequency components 
up to 300 kHz with an input of up to 1 volt (rms). The 
maximum output possible was 10 volts (peak- to-peak) with an 
expected distortion less than 0.01%. The bandpass was set 
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at DC/300 kHz throughout the experiment, thereby rendering 
the feature ineffective in comparison to the Krohn-Hite 
3550 (R) filters. An oscilloscope overlay of the input and 
reduced output of these amplifiers showed all frequency 
components to be equally amplified, with no visible 
distortion. 

When determining depth, frame reflection, and electro- 
magnetic field radiation, the transmitted pulse and received 
signal were displayed on a Tektronix 545B Oscilloscope, 
triggered by the transmitted pulse, to determine elapsed 
time of travel or amplitude of the respective signals. 

The velocimeter, described in Section III-G, sent a 
frequency signal through 150 feet of coaxial cable to a 
NUS 3010-101 Signal Conditioner and an NUS 7450-101 F-V 
Converter with an associated power supply. These provided 
frequency sampling terminals for sound velocity, temperature, 
and depth. These frequencies were converted into the values 
desired by substituting them into calibration formulas pro- 
vided. The sound velocity obtained by these formulas was 
accurate to within ±0.15 meters/sec, the temperature to 
within ±0.1° C. 

If signal processing was to be done by digital methods, 
the filtered and amplified signals were next sampled by a 
pair of Phoenix ADC 712 analog to digital converters with 
a sampling rate up to 350,000 samples per second. The 
12 bit quantization contributes an error of 1/4095 times 
the maximum voltage sampled within each frequency component. 
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The least power is generated in the 5000 Hz frequency 
component, with an average peak voltage of .050 volts. The 
error is then 1.22 x 10 ^ volts which yields a maximum 
attenuation error between the two hydrophones of 

20 L °8l0tt V 5kHz * V error )/( : V SkHz ' V error51 or - 0042 dB - 

Since the attenuation is calculated by subtracting a 

calibration run from each run in an inhomogeneous medium, 

the maximum possible quantization error is doubled, or 

.0084 dB. Because this error is random, and since 2048 

samples were taken each run, the probable error is much 

less than this value. The phase angle is also in error 

because of this quantization. The same frequency component 

yields a real and an imaginary voltage component from 

which a phase angle is computed. The error expected from 

these values is Tan ^[(V- + V )/(V , - V ___„)] 

L ^ imag error 7 v real error 

- Tan 1 [V /V ] which is a maximum when V- = 0, 
imag rear imag 

giving a value of .014 degrees. 

The digitized signal was processed finally through an 

Interdata 70 Digital Computer which analyzed the received 

signals as described in Section V-A. 

If signal processing was to be done by analog methods, 

the filtered and amplified signals were fed, separately, 

to a Hewlett-Packard 400D Vacuum-tube voltmeter and the 

individual voltage levels were recorded and compared. The 

HP 400D is an accurate average-responding rms voltmeter 

which will measure AC voltages from 0.001 volts full scale 
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to 300 volts full scale over a frequency range from 10 Hz 
to 4 MHz. The signals were then sent to a Dranetz Model 305 
Phasemeter manufactured by Engineering Laboratories, 
Incorporated. This device displayed the phase differential 
between signals received by the two hydrophones. The 
Model 305 Phasemeter operates with input voltages from 
10 millivolts to 50 volts (rms) and is accurate to ±0.25°. 

Equipment set-ups for the various experimental con- 
figurations are shown schematically in Figures 7, 8, and 9. 
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FIG. 7 EQUIPMENT SET-UP FOR DEPTH DETERMINATION 
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FIG. 8 EQUIPMENT SET-UP FOR DIGITAL SAMPLING 
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FIG. 9 EQUIPMENT SET-UP FOR ANALOG SAMPLING 



V. SIGNAL PROCESSING 



A. DIGITAL PROCESSING 

Previous research along these lines has been limited by 
the time requirements to change frequencies, change filters, 
and take amplitude and phase angle measurements at each 
frequency of interest. This required hours for a single 
series of readings before investigating other depths. The 
use of a high speed Interdata 70 Digital Computer permitted 
this series to be sampled and displayed in less than three 
minutes using the technique of the Fast Fourier Transform 
(FFT) . 

The 500 Hz sawtooth waveform sent to the transducer 
can be considered to be composed of the sum of all sinu- 
soidal harmonics of 500 Hz. The amplitudes decrease as 
the harmonics increase and there is a fixed phase angle 
relationship for each harmonic relative to the fundamental 
frequency. For example, the Fourier series representation 
of a triangular wave with a 4:1 asymmetry is 

32 °° 1 Ntt Ntt t 

f (x) = 7 Sin Sin (27) 

3tt 2 i N 2 3 L 

where 2L is the wavelength and t is the time along the 
repetitive waveform. 

A piston transducer docs not transmit all frequencies 
with equal power response or with the same acoustic 
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radiation pattern. Therefore, the signal received by the 
hydrophones is already considerably changed in form from a 
sawtooth wave, with low frequencies suppressed, and high 
frequencies unevenly attenuated by the water. The response 
of the LC-10 hydrophones further suppresses the high fre- 
quencies, and finally, the Krohn-Hite filters remove all 
frequencies below 1 kHz and above 100 kHz. Figure 2, 
page 32, shows the transmitted signal and filtered 
responses. 

This filtered, sampled signal is sent into the computer 
where a Fast Fourier Transform is performed. The FFT is 
a digital algorithm which converts the sampled waveform 
into the power contained within each frequency component 
and, unlike analog wave analyzers, preserves the phase 
relation between each frequency component and the funda- 
mental harmonic. The computer displays a power spectral 
density and associated phase angle for each of the two 
hydrophones and computes and displays attenuation and phase 
angle difference between each hydrophone. Appendix A 
shows a sample computer output generated from a contrived 
experimental run. When these results for a bubbly medium 
are compared with those for a bubble- free medium as 
expected at the 50 foot calibration depth, the change in 
sound speed and excess attenuation caused by resonant 
bubbles can be caluclated as described in Sections I-C 
and I-D. 
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The two most common problems encountered with the FFT 
are aliasing and windowing effects. Aliasing is the proper- 
ty of the FFT whereby power in harmonics with a frequency 
higher than one-half the frequency of sampling (Nyquist 
sampling rate) is combined with frequencies below this 
maximum displayed Nyquist frequency. To avoid this, we 
sampled the waveform at 320 kHz which permitted a display 
of frequencies up to 160 kHz, and we filtered out frequencies 
above 100 kHz. The windowing effect is caused by the require- 
ment for a Fourier Transform that the signal be sampled over 
a time extending from minus infinity to plus infinity. The 
FFT can be considered to take the sampled waveform in its 
finite time "window," mathematically project it into all 
preceding and subsequent time windows over all time, and 
then compute its Fourier Transform. Any discontinuities or 
radical changes in direction at the junction of these 
"windows" causes a mathematically generated series of false 
harmonics. This results in the convolution of the waveform 
sampled with a pulse having the duration of that window, or, 
in the frequency domain, this is the multiplication of 
each harmonic impulse spike with a sine function. This 
pitfall was avoided by sampling the signal at a harmonic 
of the generated signal thereby avoiding discontinuities 
and obviating the need for a hanning or hamming function. 

The 1024 sampling points permitted the FFT to compute 
512 power spectral density lines, equally divided into 
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the 160 kHz maximum frequency permitted by the Nyquist 
sampling theorem. Each power spectral density line there- 
fore represented a bandwidth of 312.5 Hz. The program 
permitted any selected frequencies to be displayed. A 
display of all spectral density lines, once with the signal 
projecting into the sound field and once with the receivers 
in the presence of ambient noise only, gave a good measure 
of the signal to noise ratio at each harmonic. Ratios in 
excess of 60 dB in the middle frequencies (35-80 kHz) 
could be obtained, and ratios greater than 40 dB existed 
throughout the spectrum. 

B. ANALOG PROCESSING 

Analog processing was used to verify the accuracy of 
the values obtained digitally, and also was used exclusively 
on those days when the computer would not operate on ship's 
power . 

Using the set-up described in Section IV-B, and shown 
in Figure 9, a monofrequency signal at each frequency of 
interest (5 kHz, 10 kHz, etc.) was transmitted through the 
F-33 transducer. A 40 volt (peak- to-peak) signal was used 
for frequencies below 80 kHz. To prevent distorting the 
transmitted signal, this was reduced to 20 volts for higher 
frequencies . 

The received signal was bandpassed to approximately 
1/10 octave, centered on the transmitted frequency, by the 
K$H 3550 filters. In order to maintain sufficient voltage 
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to drive the phase meter, signals below 25 kHz were 
amplified by a factor of 20, those "between 25 kHz and 35 kHz 
by a factor of 10, and no amplification was used for fre- 
quencies above 35 kHz. 

The amplitude of each signal was read from the HP 
400D vacuum tube voltmeter in decibels re 1 milliwatt at 
600 ohms, and the attenuation was obtained by subtracting 
the two measurements. Bubble populations were then calcu- 
lated as described in Section I-C. 

The phase difference was read directly from the 
Dranetz 305 phase meter, permitting sound speed calculations 
as described in Section I-D. 
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VI. EXPERIMENTAL RESULTS 



A. RESULTS AT ANCHOR IN MONTEREY BAY 

The first successful at-sea experiment was run on 
12 and 13 November 1974 in Monterey Bay, Lat. 36-37. 7N, 
Long. 121-54. 1W, in 38 meters of water, well removed from 
any sewage outflow or any noticeable singularities. Calm 
seas with swell height averaging 0.5 meters and wind speed 
averaging 3.5 knots, prevailed throughout the 18 hours of 
the experiment. Specific weather and sea conditions were 
recorded hourly to permit comparison with observed fluctua- 
tions. These data are tabulated in Appendix A. In-situ 
sound speed, temperature, and flow rates for each run are 
tabulated in Appendix B. 

1 . Procedure 

A 5 kHz sawtooth waveform of 30 volts peak-to-peak 
was applied to the F-33 transducer mounted on model I rack. 
Because the source has an output that is more efficient at 
mid-range frequencies, the power spectral density was not 
uniform across the frequency range of interest. Figure 10 
shows the sound pressure levels received at the closest 
hydrophone for each harmonic component of the generated 
signal. The received signals were amplified by 30 dB, 
bandpass filtered between 1 kHz and 100 kHz, then amplified 
an additional 20 dB, giving a maximum peak-to-peak signal 
of about five volts. See Figure 8 for a schematic drawing 
of the equipment set-up. 
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SPECTRAL OUTPUT OF F-33TRANSDUCER WITH 5 KHZ 
APPLIED SAWTOOTH WAVEFORM OF 30 VOLTS (PEAK-TO-PEAK) 
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FREQUENCY (KHZ) 



Two types of runs were made concurrently. One 
was an 18-hour series of runs at a constant depth of five 
meters, taken every other hour. The other type included 
two series of runs at 3, 5, 8, and 10 meters, and a refer- 
ence run at 15.2 meters, taken at 1630, 12 November and 
0800, 13 November. At each depth, five or more sets of 
samples were taken in a total time of about 15 minutes. 
These were later used to provide an average value and a 
standard deviation. Representative runs from the 18-hour 
series are included in this analysis as are each of the 
two variable depth series. 

2 . Amplitude Attenuation 

Attenuation data obtained during the 12-13 November 
at-sea experiment disproved an initial assumption important 
to our experimental analysis. The 15.2 meter depth, pre- 
viously determined to be relatively bubble- free in research 
conducted at NUC , San Diego, and Bass Strait, Australia, 
was found to have significant bubble populations at this 
location. Data analysis indicated that attenuations at 
this depth often exceeded attenuations at near-surface 
depths. Thus, there was no bubble-free reference from 
which to find excess attenuation at upper depths. It 
was decided to extract from the data the minimum dB 
difference recorded at each harmonic frequency regardless 
of time or depth and to use these values as the dB 
difference between the two hydrophones for a zero 
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attenuation medium. The computed bubble populations, then, 
are not referenced to an absolute zero bubble population but 
to a nominally bubble-free medium. 

Figures 11, 12, and 13 show graphically the number 
of bubbles in a one micron band calculated to exist in a 
cubic meter of water during this at-sea run. The shape 
and scope of these figures is remarkably similar. Popula- 
tions of small radius bubbles below 50 microns varied as R ^ 
Larger sized bubbles above 50 microns varied as R 
Some vertical stratification was apparent with more bubbles 
at shallower depths. However, this stratification trend 
was slight and did not consistently hold for all radii. 
Figures 11 and 12 depict the 1630 12 November and 0800 
13 November variable -depth series. Figure 13, the 18-hour 
survey at five meters depth, also shows only slight change 
over the observation period. However, there is slight evi- 
dence of more bubbles in the hours of darkness than daylight 

Figures 14, 15, and 16 show both excess attenuation 
and the fractional volume of air for a one micron radial 
increment contained in resonant bubbles plotted as a func- 
tion of frequency. In all three figures, a strong peak 
occurs at 15 kHz indicating a marked attenuation caused 
by bubbles resonant at this frequency. This frequency 
maximum translates to a bubble radius between 250-350 
microns. Later peaks appear to vary between 30 kHz and 
40 kHz, corresponding to bubble radii between 90 and 170 
microns . 
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Equivalent Bubble Density, n(R)dR Number/M^ dR=l micron 



FIG. 11 ATTENUATION-DERIVED BUBBLE POPULATIONS AS A 
FUNCTION OF BUBBLE RADIUS IN MONTEREY BAY 
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Equivalent Bubble Density, n(R)dR Number /M dR=l micron 



FIG. 12 ATTENUATION-DERIVED BUBBLE POPULATIONS AS A 
FUNCTION OF BUBBLE RADIUS IN MONTEREY BAY 
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Equivalent Bubble Density, n(R)dR Number /M dR«l micron 



FIG. 13 ATTENUATION-DERIVED BUBBLE POPULATIONS AS A 
FUNCTION OF BUBBLE RADIUS IN MONTEREY BAY 
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FIG 14 EXCESS ATTENUATION AND FRACTIONAL BUBBLE-TO-WATER 
VOLUME AS A FUNCTION OF FREQUENCY IN MONTEREY BAY 
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Frequency (Khz) 



FIC. 15 EXCESS ATTENUATION AND FRACTIONAL BUBBLE-TO-WATER 
VOLUME AS A FUNCTION OF FREQUENCY IN MONTEREY BAY 
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Frequency (Khz) 



FIG. 16 EXCESS ATTENUATION AND FRACTIONAL BUBBLE-TO-WATER 
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Bubble Radius, Microns 



The slight vertical stratification of bubble 
populations observed during this at-sea experiment is 
attributed to the 5 knot wind and minimal sea conditions. 

It is thought that the lack of sea surface agitation 
significantly reduced the magnitude of expected populations 
of surface origin in the upper depths. 

The presence of bubbles at 15.2 meters is attributed 
to the relatively shallow water, the in-shore conditions, 
and the probable biological activity at this location in 
Monterey Bay. 

Finally, although lessened in magnitude due to the 
use of the nominal reference, Figures 11 and 12 do show 
trends similar to those found earlier in Mission Bay 
[Medwin, 1970] as regards rate of increase of bubble popu- 
lations with decreasing bubble radius. The rate of increase 
of bubble populations with decreasing bubble radii for 

- 4 

bubbles of radii less than 50 microns agreed with the R 

found in the previous work. The number per unit volume 

- 2 4 

of larger bubbles appeared to vary as R ' compared with 
- 2 

the R rate attributable to buoyancy forces found in 
previous research. Again, there is virtually no bubble 
range showing Glotov's predicted density proportional to 
R, due to gas diffusion for R less than 60 microns. 

The deviation from trends recorded during previous 
research was initially thought to be due to test assembly 
movement in addition to the finite sampling window used in 
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digital sampling. However, the standard deviation in 
groups of five runs taken at the same depth over short 
time periods was consistently 0.1 dB or less. 

It is possible to compare the measured speed of 
sound at low frequencies with the value calculated from 
bubble populations by using equation (22). The 0200, 

13 November 5m depth plot of fractional bubble- to-water 
volume as a function of bubble radius was used. The area 
under the curve of U (R) plotted as a function of R was 
integrated over all bubble radii to arrive at a total 
fractional air- to-water volume, U. 

Equation 21 f or w much less than gives the 
asymptotic value of the differential speed 



c - c. 



= Ac = -c. 



3U 

2 (kR) 2 



In the present example, 



Ac = (-1497.18) (3/2) (1642. 5xl0 _1 °) (l/(l/36xl0 -2 ) 2 ) 



= -1.99 m/sec. 



The bubble- induced change of speed of sound at 
10 kHz for the 0200, 13 November run referenced to the 0800 
13 November 15.2m run, yields a change in sound speed of 
-1.2 m/sec. Since different references were used to 
calculate the attenuation-derived fractional volume and 
the speed of sound change as obtained from the differential 
phase, this order-of -magnitude agreement of the two values 



is reassuring. 
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3 . Sound Speed Dispersion 



During each variable -depth series a reference phase 
angle for each frequency was established using the values 
from a 50 foot run, previously considered to be bubble-free 
A three mHz velocimeter was sampled at each depth to obtain 
an accurate speed of sound which would account for changes 
in temperature, salinity, and pressure. By subtracting 
velocimeter speeds at each depth from those speeds calcu- 
lated from phase angle computations, changes in the speed 
of sound which were independent of temperature, salinity, 
and pressure were obtained. A subsequent plot of the data 
showed a requirement for an additional correction to compen 
sate for stretching of the monofilament nylon line which 
was used to pull the two hydrophones together for damping 
and for constant path length. 

The first variable -depth series was taken at 1630 
on the afternoon of 12 November 1974, in an area of 
Monterey Bay sheltered by the Monterey Peninsula from 
direct ocean swells, Lat. 36°3 7.7 N, Long. 121°54.1 . 

Wind speed was two knots, and the swell was negligible. 
Appendix B itemizes in- situ conditions for each depth. 

The differential sound speeds, after compensating for 
velocimeter values and stretch, arc shown in Figure 17. 
Figure 18 shows the standard deviation found at each 
frequency for the five runs taken at the five meter- 
depth, which proved to be typical , and for the five runs 
taken at the eight-meter depth during a later experiment 
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SOUND SPEED DISPERSION AS A FUNCTION OF FREQUENCY 
AT FOUR DEPTHS AT ANCHOR IN MONTEREY BAY 



o 

M 




O 

C\ 



o 

00 



o 



o 

vO 



o 

in 



o 

<■ 



o 

cn 



o 

e4 



74 



FREQUENCY (KHZ) 



STANDARD DEVIATION OF SOUND SPEED IN TWO DIFFERENT SEA CONDITIONS 
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FREQUENCY (KHZ) 



which had the largest standard deviation of all runs 
taken. 

The oscillating differential speeds shown in 
Figure 17 between 10 kHz and 40 kHz suggest the possibility 
of a standing waveform being transmitted by the rack at 
harmonics of the five kHz fundamental frequency. The 
intensity of this waveform, which alternately increased 
and decreased the phase angle of succeeding harmonics, 
displayed a definite depth dependence which overrode 
meaningful differential speed changes below 45 kHz. Above 
this frequency the amplitude of this effect decreased 
sufficiently to permit some observations of a change in 
sound speed with depth. At 45 kHz the speed increased 
0.35 m/sec as the depth changes from 15.2 meters to 3 meters, 
with consistently increasing values of the sampled speed at 
intermediate depths. A similar but smaller consistent 
increase holds between 70 kHz and 75 kHz. At 55 kHz there 
is a slight but consistent decrease of sound speed between 
10 meters and 3 meters. The change of direction of sound 
speed change with decreasing depth between 55 kHz and 70 
kHz, as indicated by a positive slope crossing the axis, 
suggests an increase in resonant bubble population as the 
surface was approached with a bubble resonant frequency 
of approximately 62.5 kHz. A similar presence of resonant 
bubble populations centered at 17.5 kHz and 27.5 kHz would 
be suggested except for the overriding harmonic data 
envelope at those frequencies. A strong resonant bubble 



population which increases with decreasing depth can be 
expected at 45 kHz where the harmonic interference is first 
overcome by a bubble - dominated response. Plots of excess 
attenuation versus frequency show resonant bubble popula- 
tions similar to those predicted by sound speed changes 
at 15 kHz, 40 kHz, and 65 kHz as predicted, plus peaks at 
50 kHz and 75 kHz. The absence of such a population at 
30 kHz possibly accounts for the rapid attenuation at 
that frequency of the overriding harmonic waveform. Dis- 
counting the data below 45 kHz, the maximum bubble- induced 
change in the speed of sound with changing depth at any 
one frequency is about 0.35 m/sec. Since Appendix B shows 
the velocimeter sound speed to have changed by 4.69 m/sec 
during this variable -depth series, it can be seen that 
sound speed effects caused by variation with depth of 
temperature, pressure, and salinity exceed those of bubbles 
under these conditions 

A second variable -depth series was taken at the 
same position commencing at 0800 on the morning of 13 
November. Wind speed for this series was seven knots, and 
the swell was 0.65 meters. Appendix B itemizes in- situ 
conditions for each depth. The differential sound speeds, 
after compensating for velocimeter values and stretch, are 
shown in Figure 19. Again, an overriding oscillating 
differential speed which changes in intensity with depth 
is present below 45 kHz. There is again evidence of 
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FIG. 19 SOUND SPEED DISPERSION AS A FUNCTION OF FREQUENCY 
AT FOUR DEPTHS AT ANCHOR IN MONTEREY BAY 
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excess resonant bubble populations which increase with a 
decrease in depth at 45 kHz and 57.5 kHz. A plot of 
excess attenuation versus frequency shows such a popula- 
tion at 50 kHz and 65 kHz, among others. The maximum 
bubble -induced differential speed above 40 kHz is 0.5 m/sec. 
The maximum change in sound speed shown in Appendix B for 
the velocimeter is 0.45 m/sec during this variable- depth 
series. This suggests that with sufficient mixing, the 
bubble effects on the speed of sound are on the same order 
of magnitude as those of temperature, salinity, and pressure, 
under these conditions. 

Twelve runs were taken at the constant depth of 

I 

five meters at this same location between the hours of 
1600 on 12 November and 1100 on 13 November 1974. The 
differential sound speeds of three of these runs are 
shown in comparison to a fourth run at the same depth in 
Figure 20, to give an indication of the magnitude of bubble 
induced sound speed changes at one position over a long 
period of time. These particular times were chosen to 
show the variations between sunset, nighttime, sunrise, and 
daytime. By making the reference run one of these runs, 
any depth dependent rack generated harmonics are eliminated. 

From the figure it becomes obvious that although 
frequency dependent variations are present, they are 
insignificant in magnitude in this survey. The maximum 
change in the speed of sound caused by bubbles is o.8 m/sec, 
whereas the maximum as indicated by the velocimeter during 
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FIG. 20 SOUND SPEED DISPERSION AS A FUNCTION OF FREQUENCY AT 
CONSTANT DEPTH IN MONTEREY BAY: 18-HOUR SURVEY 
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these same four runs is 3.94 m/sec. Since depth, and 
therefore pressure, was constant during these runs, this 
suggests that temperature and salinity changes near the 
surface over a 24-hour period are much more important 
factors than bubble populations under these conditions. 

B. OPEN-OCEAN RESULTS NEAR MONTEREY BAY 

A second successful at-sea experiment was conducted on 
25 November 1974. This experiment included two variable - 
depth series in areas unsheltered from ocean swells near 
Monterey Bay. The first series was taken over Monterey 
Canyon in 1600 meters of water with the vessel adrift during 
the experiment. The initial position was Lat. 36 43.5 N, 
Long. 122 05.1 W, and the position at the end of the series 
was Lat. 36 43.1 N, Long. 122 02.2 W. During the series 
the average swell height was three meters, and the average 
winds were 13 knots. The second series was taken at anchor 
near the tip of the Monterey Peninsula in 46 meters of 
water, Lat. 36 38.4 N, Long. 122 57.4 W. During this 
series, the average swell height was two meters, and the 
average wind speed was 6 knots. Specific weather and 
sea conditions are tabulated in Appendix A, and in-situ 
conditions are tabulated in Appendix B. 

1 . Procedure 

The same equipment set-up was used as described for 
the 12-13 November experiment, except that the underwater 
rack Model I was stuffed with steel wool to reduce rack 
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generated harmonic signals. An initial run at 22.9 meters 
was used for each variable -depth series to try to find a 
more nearly bubble-free reference. 

2. Amplitude Attenuation 

Probing of the 15.2 meter and 22.9 meter depths in 
search of a bubble-free reference again showed the presence 
of bubbles at these greatest depths. Even in the deep water 
run made over the Monterey Canyon, these two depths evi- 
denced bubble populations often in excess of populations 
found in the near-surface regions. The lowest attenuation 
was again used as the nominal bubble-free reference as in . 
the previous at-sea run. 

Figures 21 and 22 show the number of bubbles per 
cubic meter in a one micron band as a function of bubble 
radius for the two vertical series. As in the previous at- 
sea series conducted 12-13 November, only slight stratifica- 
tion of bubbles with depth was noted. These were unexpected 
results in view of the deeper water and the increased wind 
and sea conditions. There was less correlation in slope 
of the plot with results obtained in previous research in 

Mission Bay [Medwin, 1970]. Populations of small bubbles 

- 4 9 

below 50 microns varied as R . Larger sized bubbles 

- 2 5 

below 50 microns varied as R ’ . Figures 21 and 22 show 

marked similarity in slope but exhibit a noticeable dif- 
ference in point spread. This indicates a more pronounced 
vertical stratification of bubbles with depth in deeper 
ocean water further removed from immediate coastal effects. 
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Equivalent Bubble Density, n(R)dR Number /M dR=l micron 



FIG. 21 ATTENUATION-DERIVED BUBBLE POPULATIONS AS A 

FUNCTION OF BUBBLE RADIUS OVER MONTEREY CANYON 



7.6 NM from shore 

Swell 3.0m Wind 13 knots 

Water depth 1000m 



2000 25 November 1974 
Attenuation-derived data 
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Equivalent Bubble Density, n(R)dR Number/M dR=l micron 



FIG. 22 ATTENUATION-DERIVED BUBBLE POPULATIONS AS A 
FUNCTION OF BUBBLE RADIUS AT ANCHOR OUTSIDE 
MONTEREY BAY 



1.6 NM from shore 
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The wind and sea conditions were also more pronounced over 
the Monterey Canyon which further supports the increased 
variability in bubble populations. 

Figures 23 and 24 show both excess attenuation and 
fractional volume of air in resonant bubbles as a function 
of frequency. Strong peaks occur between 15 kHz and 20 kHz 
and again at 45 kHz. A small maximum is also evidenced 
at 75 kHz. These attenuation maxima translate to signifi- 
cant bubble populations of radii 200 to 300 microns, 80 
to 100 microns, and 50 to 60 microns, depending upon depth. 
These bubble maxima are supported by the phase analysis 
calculations used to compute sound speed dispersion. 

3 . Sound Speed Dispersion 

The differential sound speeds obtained during the 
variable depth series taken over Monterey Canyon are shown 
in Figure 25. The positive slopes crossing the axis in the 
vicinity of 40 kHz and 75 kHz indicate increasing bubble 
populations with decreasing depth. The eight-meter run 
shows a large resonant bubble population at 18 kHz, not 
seen at the other depths. A plot of excess attenuation 
versus frequency confirms that such populations exist at 
20 kHz, 45 kHz, and to a very small extent at 75 kHz. The 
greatest change in the differential sound speed at any one 
frequency was 3.01 m/sec at 10 kHz, which was obtained 
during the eight -meter run. This run consistently had a 
much larger standard deviation than normal ( see Figure 18) , 
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FIG. 23 EXCESS ATTENUATION AND FRACTIONAL BUBBLE-TO-WATER 

VOLUME AS A FUNCTION OF FREQUENCY OVER MONTEREY CANYON 
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FIG. 24 EXCESS ATTENUATION AND FRACTIONAL BUBBLE-TO-WATER 
VOLUME AS A FUNCTION OF FREQUENCY AT ANCHOR 
OUTSIDE MONTEREY BAY 
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FREQUENCY (KHZ) 



reaching a maximum of 3.66 m/sec at 10 kHz. If this run is 
discounted because of its excessive variation, the maximum 
differential sound speed is 1.4 m/sec. Appendix B lists 
the velocimeter-obtained variation of sound speed with 
depth as 0.3 m/sec, suggesting that bubble- induced sound 
speed variations with frequency under these mixed water 
conditions are more significant than those caused by 
temperature, pressure, and salinity changes. 

The differential sound speeds obtained while at 
anchor off the tip of Monterey Peninsula are shown in 
Figure 26. Positive slopes cross the axis at 17.5 kHz, 

40 kHz, and 75 kHz. Excess attenuation plots confirm large 
resonant bubble populations at 15 kHz and 45 kHz, but show 
consistently increasing bubble populations with increasing 
frequency above 70 kHz, with no unusually large populations 
around 75 kHz. The maximum differential sound speed which 
can be attributed to bubbles is 0.8 m/sec. 
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FIG. 26 SOUND SPEED DISPERSION AS A FUNCTION OF FREQUENCY AT FOUR DEPTHS 
AT ANCHOR OUTSIDE MONTEREY BAY 
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VII. SUMMARY AND CONCLUSIONS 



A portable transducer-hydrophone system with a fixed 
path length between two hydrophones of 0.95 meters was 
used to study the in-situ speed of sound and sound pressure 
amplitude attenuation as a function of frequency over the 
range of 5 kHz to 95 kHz. These experiments permitted the 
calculation of microbubble density variations with depth 
and also at a constant depth over an 18-hour sampling period 
at three different locations in and near Monterey Bay. This 
study was of interest in determining the error caused by 
omitting frequency dependence in traditional ray tracing 
methods of determining sound paths. 

The experiments were conducted in water depths of 
38 meters, 46 meters, and 1000 meters, in sea states one 
and two. Phase and amplitude measurements were taken at 
depths of 22.9 meters, 15.2 meters, 10 meters, 8 meters, 

5 meters, and 3 meters. An 18-hour survey was also taken 
at a constant depth of 5 meters in 38 meters of water. 

Analysis of amplitude attenuation in the upper depths 
of the ocean shows the definite presence and effect of 
resonant microbubble populations. However, the magnitude 
of the bubble population and its distribution throughout 
the medium is a function of many parameters. At all 
locations probed during this research work, a complicated 
dependence of bubble densities on depth was found. In some 
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cases there was greater attenuation at 15.2 meters than 
at 3 meters. Apparent peaks in bubble population, as 
evidenced by a large excess attenuation, were noted but 
showed no obvious correlation with wind or sea conditions. 
The lack of significant increase in bubble populations with 
increase in wind and sea conditions may be attributable 
to the principal attenuation sources being biological 
rather than physical. It is expected that under conditions 
of increased agitation, the near-surface bubble populations 
would be significantly higher. 

Changes in the speed of sound which could be attributed 
to bubbles were obtained by subtracting the speed of sound 
as indicated by a 3 MHz velocimeter from that obtained from 
in-situ phase angle measurements. These changes were 
plotted as a differential speed of sound for each variable - 
depth series, referenced to the speed at the deepest posi- 
tion sampled. They were also plotted for four of the 
constant -depth runs of the 18-hour survey referenced to 
one of those four runs. 

Analysis of these data indicated that the relative 
importance to the speed of sound in water of resonant 
bubbles compared to the traditional parameters of tempera- 
ture, depth, and salinity is very dependent on the degree 
of mixing of the water. There was little indication of 
any relationship between bubble- induced sound speed changes 
and wave height when all frequencies were considered. A 
larger differential sound speed occurred in the heaviest 
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sea state at the lower frequencies, with a maximum value of 
3.01 m/sec at 10 kHz. During the same series the maximum 
sound speed change caused by other parameters was only 
0.3 m/sec. 

The changes of sound speed with depth at various 
frequencies were of sufficient magnitude to suggest the 
need for a thorough in-situ study of the sound speed dis- 
persion with depth in all possible sea conditions at the 
frequencies of any proposed or existing surface ship sonar. 
Such a study would determine if the need exists for modifi- 
cation of predicted performance characteristics in surface 
duct propagation. 
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IN- SITU VELOCIMETER SPEED, WATER TEMPERATURE, AND FLOW RATES 
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Variable 25/1852 1495.18 16.75 ± 1.0 1.2 ± .4 22.9 Drifting 
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TRUNCATED SAMPLING PROGRAM OUTPUT 



H 

tD 

CX 

H 

PD 

O 






w 

H 

PD 



s 



o 



Q 

X 

w 

I! 

O 

O 

W (Vj 





O 


<>• 










<>• 


rH 




Q 












cy 


W 










II 


w 


c*: 






✓ — \ 






ps 


h-H 




<>• 


s 


Ph 


CO 


Uh 


CO 


C^* 




E 


o 


H 




w 




H 


X 




cx 


W 


Q 


H 


PD 


X 


w 




X 




PD 


Oh 


Q 


x> 


X 


Oh CO 


H 


O 


H 


S c "’ 


x 


w 


S o Q 


CO 


H 


X) 


s 


< 


XI 


< o CO 


h-H 




O 


§ VO w 


> 


Oh 


CO • Oh 


H 


i — i 




rH E 




pp 


O 




& 


Q 


rH ^ 


H 


O 


H o E 


H 


Oh 


Pi 


C-« rH H 


PD 


u 


DOU 


PO 




O 


LO H 


Oh 




Oh O HH 


Oh 


Q 


u 




X 


o o 


2NXO 




CO 


w 


H rH O 


HH 


rH ^ 


H to H 


KH 


O Oh 


rH Pi 


O Q rH U 



<N* 



97 



CONTINUE 



CHANNEL NO. 
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TRUNCATED AVERAGING PROGRAM OUTPUT 
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COMPUTER PROGRAMS 
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DIMENSION DATA (1) 
IPO=2 
IP3=IP0*N 
I 3REV=1 
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IARY (1,1 
IARY (2,1 
CONTINUE 



REQUEST INPUT PARAMETERS 
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TRANSFER TO COMPLEX ARRAY 
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OUT (1 , KK) = ARY(1,KK) 
OUT (2 ,KK) = ARY (2 ,KK) 
ARY (1, K) = TEMPI 
ARY (2 , K) = 0.0 
ARY Cl, K- ID = TEMP 2 
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501 ANGLE = 2 
GO TO 530 

502 ANGLE=0 . 0 
GO TO 530 



503 ANGLE=90 . 0 

GO TO 530 

500 ANGLE- ATAN (ARY (2,1) /ARY (1,1)) 

ANGLE = ANGLE*57. 29577951 
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WRITE (3,460) 

460 FORMAT (/13X,4HFREQ,3X,16H10LOG(PSD1/PSD2) ,3X , 13HANGLE1-ANGLE2/) 
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FIND AND READ IN SPECIFIED DTG’S 
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************************************* 

CALCULATE MEAN AND STANDARD DEVIATION 

************************************* 



Q 








H 








CO 








/— \ 




O 




Q 




» 




CO 




o 




Pu 




VO 




N— ' 




to 




CD 




+ 




O 




t — \ 




i-C 




eg 




DC 




hH 




eg 








rH 




hH 




* 




v— ^ 




X 




w 




to 




hD 








CD 




a: 




a: 




< 




< 




pc 




n 




S 




/— \ 








eg 




f — “\ 




hH 




Q /" N 








CO Q 




hH 




CP H 








^CO 




PC 




CD 




hC 




o pc 




CD 




-c i-c 




a: 




DC CD 




< 


/ ■> f N 


to 






< < PH 


rH < 




/ — 


CD PP CP 


-DC 




o eg 


> Pu, pp pp 


X cn 




/— s * hH 


< CD hH hH 


CnJ * 




eg o * 


i > Q Q 


-X 




hH • hH 


/—n<; * * 


C/vo 




#sH w 


eg i < CP 


pc - 




H |C PC 


hH /-> PP PP 


VC ^ 




^ • kC 


* eg pp pp 


pp < 




Q /->CD > > 


hH hH hH hH 


DC W 


/— \ Pc 


CO eg ^ hH hH 


w -Q Q 


^ 


*~D 


PL, hH < Q CD 


PC ph ^ /V -' 


oo * 


* 


h-j + * + W 


PD rC ^ + + 


o X PC 


H rH 


*CP hH < CP < 


-CD Q < CP 


cn lo j 


< II o 


OHS '2 S ^ O 




*\CD 


O eg • 


• II O W D CD CD * 


• II < CP DD CD 


to 


CHO 


0 

I 

s 

L 

s 

E 

s 

s 

0 


o hH CO CO PC 


w < 


PC 


CD DD 


II II DD 


H X 


O II 


II O II II ^ II II II 


ii o ii ii 


PC < o 


II rH 


cn <; i— i 


O < CP PH 


H ^ t-i 


eg <C 


cp i — i cp v — / <C £-■* CP <£, <C 


Cp eg pp pp cp E—* 


i— • ~ 


> IS 


s ^ cc x cd cd s 


2 


PC o X 


hOD)DOD)P.DO>>OOOhhDDO 


^ PC LO 


Q CD CO 


COQCOhHCOU<<COCOQQQCOCOU 


rH 









o oo 


o 


O 


e- o 


cn 


O 


rH CH 


rH 


eg 



116 



STDA = SQRT (SUMA/DIV) 

STDP = SQRT (SUMP/DIV) 

WRITE (3,910) FREQ (1, Iw) ,AVGP, STDP, AVGA, STDA 
910 FORMAT (IX ,F8 . 1 ,4 (1X,F14. 5) ) 

210 CONTINUE 
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